Thermodynamics of 0{N) sigma models: corrections 
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The thermodynamics of the 0{N) linear and nonlinear sigma models in 3+1 dimensions is studied. 
We calculate the pressure to next-to- leading order in the 1/A'' expansion and show that at this order, 
temperature- independent renormalization is only possible at the minimum of the effective potential. 
The expansion is found to be a good expansion for A'^ as low as 4, which is the case relevant for 
low-energy QCD phenomenology. We consider the cases with and without explicit symmetry break- 
ing. We show that previous next-to-leading order calculations of the pressure are either breaking 
down in the temperatures of interest, or based on unjustifiable high-energy approximations. 



I. INTRODUCTION 



It is well-known that although the QCD Lagrangian 
possesses a chiral symmetry in the limit of zero quark 
masses, the true QCD ground state does not respect this 
symmetry. The chiral symmetry is spontaneously broken 
by quantum effects. To be specific, the QCD Lagrangian 
with Nf massless quarks has a global SU{Nf)L x 
SU{Nf)}i symmetry, which for the ground state at low 
temperatures is broken down to an SU{Nf)v symmetry. 
According to Goldstone's theorem, there is a massless, 
spinless particle for each generator of a broken global 
continuous symmetry. In this case this implies the oc- 
currence of Nj — 1 Goldstone bosons. In phenomenolog- 
ical applications Nf is either two or three, and one also 
has to take into account the explicit symmetry breaking 
due to the nonzero quark masses. Both the spontaneous 
and the explicit chiral symmetry breaking are apparent in 
the low-energy hadronic particle spectrum, where the ex- 
pected number of relatively light mesons is observed (e.g. 
the three pions for Nf = 2). At sufficiently high temper- 
atures one expects the chiral symmetry to be restored 
and lattice simulations of QCD suggest that this happens 
at a temperature of approximately 150 MeV depending 
on the number of quarks and their masses. Heavy-ion 
collisions at RHIC and LHC are expected to reach such 
temperatures and will allow experimental studies of the 
deconfined, chirally symmetric phase of QCD Q . 

Apart from using lattice simulations it has not yet 
been possible to calculate thermodynamic properties, 
such as the pressure, from QCD in the low-temperature 
hadronic phase. However, this can be done using low- 
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energy effective theories. Such effective theories for 
the low-energy particle spectrum, involving both mesons 
and baryons and displaying the above-mentioned pat- 
tern of chiral symmetry breaking, were constructed be- 
fore QCD itself. In the case of two flavors, the situa- 
tion is the simplest, since one can exploit the fact that 
the SU{2)l X SU{2)ji symmetry is locally isomorphic to 
SO{4:). If baryons (nucleons in this case) are not in- 
cluded the simpler 0(4) linear sigma model can be used 
as a low-energy effective theory for describing the dy- 
namics of three pion flelds and one sigma field. These 
four fields form a four-dimensional vector (f) in the fun- 
damental representation of 0(4). At low temperature, 
the 0(4) symmetry is spontaneously broken down to 
0(3), where the sigma field acquires a vacuum expec- 
tation value and the three pions are interpreted as the 
Goldstone bosons. For Nf > 2 there is no connection 
between the SU{Nf)L x SU{Nf)R model and the 0{N) 
linear sigma model, but the latter has been studied in 
great detail for general N due to its relevance to spin 
models. 

In this paper the 0{N) linear sigma model (LSM) 
and 0{N) nonlinear sigma model (NLSM) in 3-1-1 di- 
mensions will be studied at finite temperature and to 
next-to-leading order in the 1/A^ expansion. At zero 
temperature, the 1/N expansion was applied to 0{N_ 
sigma models a long time ago at leading order (LO) ^ 
and at next-to-leading order (NLO) Q. At finite tem- 
perature the LO in 1/N has been studied in Ref. 3- 
The effective potential in the large- A/" limit is that of an 
ideal gas and thus straightforward to compute. The NLO 
1 /N corrections to the free energy involve a momentum- 
dependent self-energy and cannot be evaluated analyti- 
cally. In Ref. 5] , the author therefore carried out a high- 
temperature expansion to obtain purely analytical results 
for the LSM. Similarly, in Ref. the authors resorted 
to a "high-energy" approximation, which makes the cal- 
culations manageable. However, this approximation is 
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uncontrolled and it is difficult to assess how reasonable 
it is, unless one calculates the full NLO corrections. 



case, since the second-order phase transition turns into a 
smooth cross-over. 



The 0{N) sigma models have also been studied 
in detail at finite temperature using various other ap- 
proaches. A systematic study has been carried out by 
Chiku and Hatsuda using optimized perturbation the- 
ory . The method was used to calculate spectral func- 
tions, properties of the effective potential, and dilep- 
ton emission rates. This method is convenient from the 
point of view of renormalization. The Cornwall- Jackiw- 
Tomboulis or 2PI formalism ^ has also been used to ex- 
amine various properties of the 0(N) linear sigma mod- 
els at finite temperature 0, 0, O) 1^ ^i ^3, , see 
Ref. for a recent review. For example, in several pa- 
pers the temperature dependence of the pion and sigma 
masses, and of the vacuum expectation value of the sigma 
field, have been investigated. In the low-temperature 
phase, the 0(N) symmetry is spontaneously broken and 
it is expected that the symmetry is restored via a second- 
order phase transition. The calculations of the effective 
potential as a function of temperature have been car- 
ried out in the Hartree approximation and the large-iV 
limit [H in El El El lli. In these cases, the gap 
equations for the propagators are easy to solve since the 
self-energy reduces to a local mass term. In the Hartree 
approximation, the result has been shown to be problem- 
atic (and a first-order phase transition occurs), which has 
been remedied by including more diagrams in the trun- 
cation .12» .18], resulting in a second-order phase transi- 
tion. If one goes beyond the Hartree approximation or 
includes the NLO contributions in the 1/A^ expansion, 
the gap equations become nonlocal and very difficult to 
solve. 



The NLSM in 3-1-1 dimensions is nonrenormalizable 
and should be viewed as an effective theory, which is 
valid up to a certain energy scale where new physics en- 
ters. Strictly speaking the LSM is renormalizable, but 
since it becomes a trivial theory in the limit where the 
cutoff goes to infinity, we treat it as a theory with a fi- 
nite cutoff. Given a finite cutoff, we speak of divergences 
when terms are increasing in magnitude without bound 
as the cutoff is increased. The low-energy physics should 
be independent of such terms (decoupling) and one can 
subtract them in the renormalization procedure in order 
to avoid increasing sensitivity to the ultraviolet cutoff as 
it grows. On general grounds, one expects the tempera- 
ture dependence to be insensitive to an increasing cutoff 
due to the exponential suppression provided by the Bose- 
Einstein distribution. Therefore one expects the renor- 
malization to be possible in a temperature-independent 
way. 

The paper is organized as follows. In Sec. II, we dis- 
cuss effective actions of the LSM and NLSM in the 1/A^ 
expansion. In Sec. Ill, we calculate the effective poten- 
tial and gap equations at NLO. In Sec. IV, we discuss our 
results for the pressure at NLO for general N and for the 
special case of = 4. Also, the so-called high-energy 
approximation is discussed and compared with exact nu- 
merical results. In Sec. V, we elaborate on the choice of 
parameters for A^ = 4, in order to make contact with low- 
energy QCD phenomenology. In Sec. VI, we summarize 
and conclude. 



In this paper, we will study the IPI effective poten- 
tial at NLO in the 1/N expansion without resorting to 
a high-temperature or high-energy approximation. Thus 
the analysis presented here is an extension of the papers 
by Jain p , and by Bochkarev and Kapusta . We will 
follow the approach to the NLSM in l-fl dimension of 
Ref. 01 . In the present case wc do not need to deal 
with thermal IR renormalons, which simplifies the renor- 
malization procedure. Nevertheless, we will obtain simi- 
lar conclusions about the renormalization of the effective 
potential in 3-t-l dimensions as in 1+1 dimensions. It 
turns out that at NLO, temperature-independent renor- 
malization is only possible at the minimum of the effec- 
tive potential. This aspect of the 1/A^ expansion was 
missed in previous work d |^ , since the renormalization 
is considerably simplified or even ignored in the various 
approximations. 

Since explicit chiral symmetry breaking plays a very 
important role in the actual hadron spectrum at low en- 
ergy, we will consider also the case of explicit symme- 
try breaking. The results change considerably and more- 
over, a critical temperature cannot be determined in that 



II. EFFECTIVE ACTIONS 



The Euclidean Lagrangian of the 0(A^)-symmetric 
linear sigma model with a symmetry breaking term H is 
given by 
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kc^,-Nfl,f -VnH4>n ,(1) 



where i — 1 . . . N. Summation over repeated indices is 
implicitly understood. The subscript b denotes a bare 
quantity. The coupling constants are rescaled with fac- 
tors of A^ in such a way that for large A^ the action nat- 
urally scales as A^. 

It is possible to eliminate the quartic interaction 
term from Eq. by introducing an auxiliary field which 
is denoted by a, in order to allow for Gaussian integra- 
tion. To this end, we add to the Lagrangian Eq. the 
term 



(2) 
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such that one has 



N o ^ 



(3) 



Using the equation of motion for a, one recovers the orig- 
inal Lagrangian in Eq. JQ). In the limit Af, — s- oo, one 
obtains the Lagrangian of the nonlinear sigma model. 

If explicit symmetry breaking is absent (H — 
0), the field (j) acquires a vacuum expectation value 
by spontaneously breaking the symmetry. Because of 
the residual 0{N — 1) symmetry, we can write (f> = 
(tti, 7r2, . . . , ttat^i, (t), such that only ~ c has a 

nonzero expectation value. For H > the same argu- 
ment applies, because the action is minimal when the a 
field is the only one that acquires an expectation value. 

Integrating over the tt's gives the following effective 
action 

Soft = i(iV-l)Trlog(-a2-m) 



dr / d^x 



2 id,<^r 



(4) 



We next parametrize the quantum fields a and a by writ- 
ing them as a sum of space-time independent vacuum 
expectation values irri^ and ct, and quantum fluctuating 
fields a and a: 
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(5) 
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Using Eq. Q one can show that the vacuum expectation 
value of a is purely imaginary. The vacuum expectation 
value of a is proportional to \/N, which follows from 
Eq. CJ. Substituting Eqs. and 1© into Eq. Q), the 
effective action ScS can be written as 
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Expanding Eq. Q in powers of l/\/N up to corrections 



of order l/y/N, one finds ^ 

Scff 1 

2 



2Ah 



NHa + \/ N X terms linear in a and a 



P'^ + m 



2 , ^2 )X\ (8) 



where — (ap, ^p) is a vector containing the Fourier 
transforms of a and a, and the function n(P, to) is given 
by 



n(P, to) ^ q2 + ^2 (p + Q)2 + ^2 

We have introduced the sum-integral 
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i3o=27rnT 
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which involves a summation over Matsubara frequencies 
go and an integral over three-momenta q. For later con- 
venience we also introduce a symbol for the difference 
between a sum-integral and an integral 
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(11) 



where / = /d'*Q/(27r) 



III. EFFECTIVE POTENTIAL AND GAP 
EQUATIONS 



One can obtain the effective potential through next- 
to-leading order in the 1/A^ expansion from Eq. © by 
performing the Gaussian integral over the fluctuating 
fields a and a. Up to corrections of order the effec- 
tive potential can be written as 



V(TO2,a) - iVVLo(m2,a)-f VNLo(m2,a) , (12) 



where 



1 



log(p2 + TO^) 



Vnlo , a) 



-i^log/(F,m) 



(13) 
(14) 



^ The terms that are linear in a and a vanish at the minimum of 
the effective potential. 
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Here, 



/(P, m) = 16TT^n{P, m) 



P2 + m 



(15) 



To derive the effective potential, we subtracted divergent 
constants which are independent of ct, m and the tem- 
perature. Equivalently, these terms can be removed by 
adding a vacuum counterterm AS to the effective poten- 
tial. In the following, we simply drop these terms. 

In thermodynamic equilibrium, the system will be 
in the state that minimizes the effective potential with 
respect to m? and maximizes it with respect to a. This 
difference is due to the fact that the vacuum expectation 
value of a is imaginary and that of a is real. Differenti- 
ating the NLO effective potential with respect to and 
a, one obtains 
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2m^ 



p2 



dn(P,m) 
1 y—f dm^ 

A^Tpn(P,m)4 
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%P^+m^ n(P,m) + ^ + p|^ 

= H 



25-2 



(16) 



(17) 



These equations are often referred to as gap equations. 
Solving the gap equations gives m and a as a function of 
the parameters H and A, and of the temperature. 

The inverse a and a propagators can be obtained 
from Eq. JSJ. One finds 



Dl\P,m) ^ n(P,m) + — + 



1 



D^\P,m) = P^ +m'^ 



Xb P2 + m? 
n(P,m) + 2/Afc 



, (18) 
• (19) 



The values for vn? and u are determined by solving the 
gap equations. Using the a propagator, one finds that 
the inverse tt propagator is given by 



D~\P,m) = P'^ + 
N% {P + Qf+m^ n(Q,m) + f + 



(20) 



From this equation and the gap equation H17|l one sees 
that also at NLO in the broken phase where a ^ Q (for 
H ^ 0), the pions are massless, in accordance with Gold- 
stone's theorem. 

From Eq. (|19|l one can see that in the unbroken 
phase, the a mass becomes equal to the LO mass of the 
TT field, which is m?. It may appear therefore that the 
a and tt masses are not equal at next-to-leading order, 
but this is not a correct conclusion. We note that the a 



field only starts to propagate at NLO, so its 1/N mass 
corrections require a next-to-next-to-leading order calcu- 
lation. We emphasize that in the calculation of the NLO 
pressure one only needs the LO masses, as we will see 
explicitly below in Eq. H53() . 



In subsecs. IIII Al and IIII Bl we will explicitly calcu- 
late the leading-order and next-to-leading order contri- 
butions to the effective potential in 3-1- 1 dimensions. We 
will evaluate integrals using an ultraviolet momentum 
cutoff A and assume that A ^ m, 2ttT. 



A. Leading-order contribution 



The leading-order contribution to the effective po- 
tential is 
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Jo (/3m) 
+ Ha .(21) 



Eq. H21|l contains ultraviolet divergences in the sense ex- 
plained in the introduction. These divergences can be 
dealt with by defining the renormalized parameters 
and A as 



/2 = fl,-A'/16n^ 
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= log 



32w' 



(22) 
(23) 



where A = A(/x). After this renormalization, the leading- 
order effective potential becomes 



647r2 



Vlo = ^{f.-^' 
where the function Jo(/9m) is 



327r2 



Jo{f3m) = 



32 

3T4 



dp — n{ujp) 

UJri 



1 

~ 2 
(24) 



(25) 



Here, n{ujp) = [exp(/3u;p) — 1]^^ is the Bose- Einstein dis- 
tribution function. Note that one makes an error in the 
evaluation of Jq by integrating up to infinite momenta in- 
stead of up to A. However, this error is negligible as long 
as A m, 2ttT. This remark also applies to the func- 
tions Ji, Kq and A'j^ defined below. For an investigation 
of how to deal with a finite cut-off in the calculation of 
sum-integrals cf. Ref. [20| . 

The renormalization group equation for the running 
coupling A that follows from Eq. (|23|l is 



/3(A) = M 



A2 



dX 



(26) 
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Note that the /5-function is exact to all orders in in 
the large- -/V limit, but differs from the perturbative one 
obtained at one loop. However, at NLO they agree. 

Since the potential term in the Lagrangian should al- 
ways have a minimum in order to have a stable theory, 
must be positive (cf . [ij for a detailed discussion) . From 
Eq. H23() it immediately follows that there is a maximal 
value for the cutoff given by 



Amax = M exp 



(27) 



Therefore this theory should be viewed as an effective 
theory, which is valid up to the cutoff given by Eq. (|27|l . 
Taking the cutoff to infinity is equivalent to taking A 
to zero, which implies that the theory is trivial. One 
should keep in mind that the renormalized leading-order 
effective potential does not depend explicitly on A, but is 
only valid for m and T much smaller than Amax- When 
A = Ajjiax the linear sigma model reduces to the nonlinear 
sigma model, since in this case Xi, = oo. 

The leading-order renormalized gap equations follow 
from differentiating Eq. 124|) with respect to and a 
and are given by 



G = WTT^f^ , 



(28) 
(29) 



where 



327r2m2 



G = rVi(/3m) + le-K^a^ - log 

Here, we have defined the function Ji(/3m) as 
8 f°° 

Ji{(3m) = 7^ / dp — n(wp) . 
Jo 

If i7 = 0, one can show by using the gap equation H29() 
that either to = or a = 0. From the gap equation H28|) 
it follows that for m = the expectation value of a has 
the temperature dependence 



(30) 



(31) 



2^2 



(32) 



Ax T = Tc = V 12/77 there is a second-order phase tran- 
sition Isj. Below Tc the 0{N) symmetry is broken spon- 
taneously to 0{N - 1) since a ^ 0. Above the 0{N) 
symmetry is restored and one has a = and to 7^ 0. 



B. Next-to- leading order contribution 



In this section we will show that it is not possible 
to renormalize the next-to-leading order effective poten- 
tial in a temperature-independent way. ft turns out that 



we can only renormalize the effective potential at the 
minimum, since the temperature-dependent divergences 
become temperature independent by using the leading 
order gap equations. To show this, we will extract the 
divergent parts of the effective potential, which can be 
done analytically. 

In order to isolate all divergences, in principle we 
need to evaluate H(P, to) including corrections of order 
TO^/A^, since such terms can also give rise to divergences 
in the effective potential. However, since the linear sigma 
model is an effective theory, Eq. |^ should be viewed 
as the part containing only the relevant operators. For 
instance, we have not included irrelevant operators of di- 
mension six, which also contribute to n(P, to) at order 
1/A^. Therefore, for consistency with Eg. we do not 
consider order 1/A^ terms in H(P, to) [21|. We only re- 
tain the unsuppressed terms of H(P, to) in an expansion 
in 1/A^, since this expression is much less complicated 
than the exact one. We find 



H(P, to) 
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log 
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167r2 





P2 



P2 



4m2 / VP^ 

log 



4to2 



/p2 



VP2 



4to2 + V P2 

-FHt(P,to), (33) 
where the temperature-dependent part of n(P, to) equals 



Ht(P,to) 
1 



q ( q' +pq + A'^\ 
— log ( — r— ^ ) n{ujq) . (34) 



Here 



Stt^P Jq ujq \q^ - pq + 



,2 _ P^ + 4m2p2 
4P2 



(35) 



In the limit P ^ m,T , we can approximate T\t{P, w) by 
I 



Ht(P,to) 
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T2 Am'^T'^r? 



(3P2 - 4p2) 

P6 



Jo(/3to) 



(36) 



The next-to-leading-order effective potential has 
only ultraviolet divergences. Using the leading order 
renormalization of Af,, it is easily seen that /(P, to) be- 
comes finite. Also, the difference 



^ log /(P, to) 



(37) 



is finite (cf. Sec. lIVBl) . Therefore, all possible divergences 
of Vnlo can be isolated by calculating 



J log/HE(-P,m-) , 



(38) 
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where /he (-P, w) is the high-energy (HE) approximation 
to I{P, to). It gives the large-P behavior of I{P, m). Af- 
ter averaging over angles, we find 



1 r 1 ^ 1 ^2 I (C2 



P4 Cl 



where 



^2 = 



1 , \ , 327r 
log ( ^ I + 1 



A 

-2to^ 1 + log 

\TO^ y 

I ™4 



21og(5)-l 



[327r2CT2 + 2TVi(/3?: 



(39) 



(40) 



(41) 



(42) 



By integrating the function log /he over P, we obtain all 
the divergences of the NLO effective potential. The loga- 
rithmic and power divergences are given by the quantity 
/), which is 
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gl+327r2/At, 

1 + 2ei+32-VA. li 
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.l+327rVA[, 



2mnog( ^ 



(43) 



while the terms that have a small cutoff dependence 
through their dependence on Afc, are given by the quan- 
tity which is 
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327r^ 
Ah 



(44) 



Since G depends explicitly on the temperature, it is im- 
possible to renormalize the next-to-leading-order effec- 
tive potential in a temperature-independent way. How- 
ever, at the minimum, one can use the leading-order gap 
equation to show that G = IGtt^/^. Hence, the 

divergences become independent of the temperature at 
the minimum and we can renormalize in a temperature- 
independent manner. We discuss this next. 

The divergence in the first line of Eq. H43I) is in- 
dependent of TO in the minimum. This divergence can 
be removed by vacuum renormalization. The divergent 
terms which are proportional to rri} can be removed by 



defining the renormalized parameter /^r as 

A2 



f2 =f2 

J TT J TT.b 
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,l+32ir2/A6 



(45) 



The remaining divergence is proportional to to^ and is 
removed by renormalizing A;, as follows 

32^2 32^2 / 8\, /A2^ 
■ ^ I 1 + — I log I — 



A 



Xb 



N 



log 



(46) 



From Eq. (|46(l . we obtain the /?- function governing the 
running of A: 

m -1^(^ + 1) ^ (47) 



167r2 



which coincides with the standard one-loop /3-function 
in perturbation theory. One could argue from the renor- 
malization that one can only trust the 1/iV expansion for 
iV ^ 8. Although the 1/A^ correction to the /3 function 
indeed has a large coefficient, this is not the case for the 
effective potential itself as we will see below. 

As mentioned, the terms in E have a small cutoff 
dependence through their dependence on A;,. We do not 
renormalize them, since they do not grow without bound 
with increasing cutoff and are not strictly speaking diver- 
gences. The effective potential does not become increas- 
ingly sensitive to them with increasing cutoff 2. The term 
in the first line of Eq. H44|) becomes smaller if we increase 
A and the absolute value of the other term from E in- 
creases as a function of A, but is bounded by a finite 
number which is independent of Ah and A. Moreover, 
renormalizing these terms would invalidate the ex- 
pansion, because of their magnitude. This is similar to 
ordinary perturbation theory, where one is only allowed 
to do finite renormalizations that do not invalidate the 
perturbative expansion. A final reason for not renormal- 
izing these terms is the connection with the nonlinear 
sigma model (Af, = oo). In that case, the terms from E 
are not divergent and we would choose to renormalize /^r 
just as in Eq. (|45|l with Xh — oo. 

The NLO correction changes the critical tempera- 
ture Tc- Since the NLO gap equations are complicated, 
we are not able to get an analytical expression for the 
NLO critical temperature. In the limit of small Xb and 
H — the gap equations however simplify to 



p2 



TO- 



2m2 
Xb 
1 



Since A < Amax, one also has to 
considering decreasing values of A. 



(48) 
(49) 

accordingly, by 
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From the gap equations it follows that the critical tem- 
perature at NLO is 



12 



1 + 2/N 



(50) 



This result is the same as obtained in Refs. This 
result is probably only correct in the weak-coupling limit 
and Tc may depend on A at NLO in 1/A^. As we will see 
in the next section the transition is of second (or higher) 
order. 



/tt = 47 MeV (note that our f^^ is 1/2 times the more 
conventional definition) and if there is explicit symmetry 
breaking H = (104 MeV)'^. Otherwise H vanishes. A re- 
alistic choice of parameters would allow us to compare to 
lattice-QCD simulations of the Nf = 2 case for T <Tc, 
although this would require extrapolation of the lattice 
results down to the actual pion masses. 



A. Leading-order contribution to the pressure 



IV. PRESSURE 



The pressure V{T) is equal to the value of the ef- 
fective potential at the minimum at temperature T mi- 
nus its value at the minimum at zero temperature. As 
we showed in the previous section, we can renormalize 
the effective potential at the minimum. The pressure is 
therefore a well-defined quantity. In order to determine 
the NLO effective potential in the minimum, we need 
the gap equation only to leading order 0. Writing the 
solutions to the gap equations as 



'"■LO 



- 0-NLo/N , 



(51) 
(52) 



and Taylor expanding the effective potential H12I) , wc ob- 
tain (up to 0{1/N) corrections) 



V(m2,CT) 



~™NLO 



+fNLO 



LO'^LO) + VNLo("tLO'^Lo) 



d'm? 

9VLo(a) 



da 



(53) 



The last two lines of Eq. H53|l vanish by using the leading- 
order gap equations. In the following, we will write the 
pressure V as 



V = NVlo + T'nlo • 
From the discussion above, it follows that 



(54) 



-Plo = VZo(m|,aT)-V^<J°(m2,ao) , (55) 
^NLo = Vnlo("^t:^t) -V^i7(?] (771^,0-0 ) , (56) 

where rri^ and ax are the solutions of the leading-order 
gap equations and if^ . at temperature T. 

In the following, we will present the results for the 
numerical evaluation of the leading and the next-to- 
leading order contributions to the pressure for general 
A''. In subsecs. IIV CI and IIVDI we specialize to A^ = 4. 
As we will motivate in Sec. we will use the follow- 
ing values for the parameters: A(/i = 100 MeV) = 30, 



Using Eqs. (|24|l and H55|l . we can easily calculate the 
leading-order pressure. In Fig. ^ we show the leading- 
order pressure normalized by T"*. If _ff = 0, the pions 
are massless below Tc and the leading-order effective po- 
tential (|24|) reduces to the free energy of an ideal gas of 
massless particles: Vlo = 



0.12 
0.1 



Vlo 

rp4 



0.08 - 
0.06 - 
0.04 - 
0.02 - 




"1 — I — I — r 



_L 



J I I L 



_ff = 

7J ^ 

I I L 



50 100 150 200 250 300 350 400 450 500 

T (MeV) 

FIG. 1: Leading-order pressure Plo normalized to as a 
function of temperature without and with explicit symmetry 
breaking. 



B. Next-to-leading order contribution to the 
pressure 



To calculate the next-to-leading order contribution 
to the pressure, we decompose T'nlo as follows 

T'nlo = Dimr) - D{mo) + + F2 , (57) 

where Dim) is the term containing logarithmic and 
power ultraviolet divergences given in Eq. (|43|l . and Fi 
and F2 are finite terms defined below. 

The term Fi has a weak cutoff dependence and is 
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defined by 




n(P,mT) + |- 



p2 



log 



n(P,mo) + — 



2 

a; 



2al 



p2 



-D{mT) + Dimo) 



(58) 



We calculated Fi numerically by rewriting the terms in- 
volving D as an integral. We can then subtract the inte- 
grands, instead of the large values of the integral. In this 
way, it is easier to avoid large numerical errors. 



The function F2 is defined by 
^ *' ' n(P, mr ) 



2 



2cr^ 



p2 



In order to calculate the function _F2j we have modified 
the Abel-Plana formula [23 to obtain the relation 

Ajv = ^ /(po = 27rr7,) ~ ^ j '^PofiPo) 

n=N ^ 
1 f°° 1 

= - / dplm/(a + *p)— — , (59) 

where a = 2-k{N — 1/2). This formula is valid as long 

/(po) has no poles or cuts for Rc(po) > a, /(Po) G 
for po G and /(po) grows slower than an exponential 
for po ~* 00. This relation is useful for numerical calcu- 
lations, since it prevents us from subtracting two large 
quantities. In our case /(po) — log /(P, to) which is even 
in Po and has a cut for Re(po) = 0. Then we can use that 



2A 



N - 



n=-N+l 



27rn)-^ I dpofipo) , 



(60) 

where > 1 because of the cut. We checked that chang- 
ing TV has no effects on the results. After calculating 
A_oo, we integrate over three-momentum p up to Amax, 
which gives a finite result for F2. We observe that the 
difference of a sum-integral and an integral is dominated 
by the low-momentum modes. This shows that the high- 
energy approximation in Ref. Q applied to their "in- 
teraction pressure" (implicitly defined as such a differ- 
ence and directly related to P2) is invalid, since the high- 
momentum modes are assumed to give the main contri- 
bution. 

After renormalization, we find that the next-to- 
leading order contribution to the pressure is 



'Pnlo 



™T 1 

8^^°S 



'"0 , 



Fi 



F2 , 
(61) 



'Pnlo -0.04 - 
-0.06 - 
-0.08 - 
-0.1 - 
-0.12 




H^Q, 

J? = 0, 

J? = 0, 



A = 10 GeV 
A = 5.0 GeV 
A = 10 GeV 
A = 5.0 GeV 



_L 



J I I I I I I L 



50 100 150 200 250 300 350 400 450 500 

T (MeV) 

FIG. 2; Next-to-leading order contribution to the pressure 
normalized to T^, as function of temperature for H — and 
H = (104 MeV)^, for different values of the cutoff A. 



which is shown in Fig. [3 At T = 0, we can calculate 
'Pjqho/T^ exactly and use it as a check of the numeri- 
cal calculations. At T = 0, clearly Pi = 0, and hence 
Vnlo/T-^ ^ F2/T^. At P = 0, it is easy to see that for 
low P, I{P,m) is dominated by 327T'^a'^ /{P'^ + m'^). This 
gives 



P2 



(62) 



For H = 0, the mass m vanishes and so Pnlo/?"* = 
-7rV90. For H = (104 MeV)^ we have that Pnlo/T'^ = 
0, since T is in that case much smaller than to, such that 
the pressure is exponentially suppressed. 

The pressure ioi H — Q is approaching the H ^ Q 
pressure at high temperatures, indicating that the ef- 
fects of the explicit symmetry-breaking terms become 
smaller at higher temperatures. This is because H is 
a temperature-independent constant. 



C. Pressure of the 0(4) linear sigma model 



In order to make contact with two-flavor low-energy 
QCD, we specialize to A^ = 4. In Fig.O we show the pres- 
sure for A^ = 4 and if = to next-to-leading order as 
function of T normalized by T^. The LO pressure below 
Tc equals the pressure of a gas of four massless noninter- 
acting scalars. This follows immediately from Eqs. H24fl 
and H25|) . At NLO the sigma field becomes massive. For 
temperatures much lower than ma-, the contribution to 
the pressure from the sigma is Boltzmann suppressed and 
we have (to good approximation) V — 7r^T'*/30, which 
is the pressure of a gas of three massless noninteract- 
ing scalars. From the calculations we conclude that the 
transition to NLO is of second (or higher) order since the 
derivative of the pressure is not diverging. 
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0.44 



"1 — I — I — I — I — I — I — I — r 




LO 

LO + NLO 

J I I I I I I I 1_ 



50 100 150 200 250 300 350 400 450 500 

T (MeV) 

FIG. 3: LO and NLO pressure normalized to T"*, for = 4 
as a function of temperature, for H = and A — 5.0 GeV. 



In Fig. ^ we show the pressure for = 4 and 
H = (104 MeV)'^ to next-to-leading order as function 
of T normalized by T'*. 




50 100 150 200 250 300 350 400 450 500 

T (MeV) 



FIG. 4: LO and NLO pressure iov N = A normalized to T"*, as 



a function of temperature for H 
GeV. 



(104 MeV)'' and A = 5.0 



In Figs. 13 and 01 we have chosen the cutoff A = 5.0 
GeV. Because we wish to make contact with low-energy 
QCD, a few comments on this choice are in order. For 
the low-energy chiral Lagrangian, the cutoff is usually 
taken to be Stt/^ (using our definition of /^), which is 
around 1.2 GeV. However, for the present purpose this 
value would be at the limit of applicability, since the crit- 
ical temperature at which chiral symmetry is (approxi- 
mately) restored is only about a factor of 8 smaller and 
we have to satisfy the requirement that 2'kT ^ A. In this 
way one ensures that one sums over sufRcient Matsubara 
modes ^ . Therefore, we have taken the cutoff consider- 
ably larger to reduce the sensitivity to the cutoff, but 



emphasize that only for the region below Tc can one ex- 
pect the result to be of relevance for the QCD pressure. 



D. Pressure of the 0(4) nonlinear sigma model 



In the limit Ah = oo, we obtain the Lagrangian for 
the nonlinear sigma model. So there are no counterterms 
for logarithmic divergences. We will only renormalize 
/7r,b as in Eq. H45() with A;, = oo. This implies that F2 
has a small cutoff dependence. 

In Fig. [5I we show the pressure of the 0(4) non- 
linear sigma model without explicit symmetry breaking 
[H = 0), through next-to-leading order in 1/iV. We have 
calculated the pressure for different values of the cutoff. 
The LO result for A = 20 GeV is included. For compari- 
son, we also show the pressure resulting from the approx- 
imations employed by Bochkarev and Kapusta, Ref. 
A considerable difference between our results and those 
of Bochkarev and Kapusta is observed. 



V_ 

rp4 



0.45 

0.4 
0.35 

0.3 
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0.15 h 

0.1 



LO, A = 20 GeV 
A = 50 GeV 
A = 20 GeV 
A = 10 GeV 
Eq. (52), Ref. [6] 
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FIG. 5: NLO pressure of the nonlinear sigma model for A'^ = 4 
normalized to T*, as a function of temperature for different 
values of the cutoff A. For comparison we have included the 
LO pressure and a curve corresponding to the NLO pressure 
expression from Ref. 0. 



We next discuss the approximations made in Ref. ^ . 
The T = part of sum-integrals are omitted such that 
every ^ is replaced by Hence H(P, m) is replaced by 
Ht(P, m). In the high-energy approximation the latter 
is approximated by 



HT(P,m) 



327r2A2 



(63) 



^ The amount of Matsubara modes to be summed o ver depends 
on how one implements the cutoff on sum-integrals l.20| . 



where Ji and are given by Eqs. H31|) and H35|l . respec- 
tively. The term involving in Eq. H15|) is omitted and 
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the pressure reduces to 



V 



--4 log- 



(64) 



2jp""'=P4 + 4TO2p2 • 

Defining the functions 

^o(P^) = ^4l dp£-n(c.±), (65) 



where uj± — \/p'^ + fn? ± m, the pressure becomes 



NT'^JoiPm) TT 



647r2 



90 



64^ 



[Xo+(/3m)+Xo-(/3m)] 



(66) 



Expanding Eq. (|66|l in powers of m/T and rescaling with 
factors of N, one obtains Eq. (52) of Ref. 

For completeness, we also give the gap equations in 
this approximation: 



167rV^ = rVi(/3m) + IGtt^ct 



2^2 



[K+{f3m)+K^{f3m)] , (67) 



m^a = 



where the functions Ki are 



8 



dp 



rp2 



n{uj±) . (69) 



There are several problems with the approach in 
Ref. 1^ . Firstly, it is incorrect to ignore zero-temperature 
contributions to the pressure and it also obscures renor- 
malization issues. 



Then, as we saw in Sec. IIVBI the 
arguments for applying the high-energy approximation 
are not valid. Furthermore the term proportional to a is 
not suppressed by 1 /N compared to the other terms ap- 
pearing in I(P,m). Fourthly, as the solutions to the gap 
equation (|67() indicate for T > Tc m/T becomes signifi- 
cantly larger than one, hence the m/T expansion breaks 
down. If one were to use Eq. (|66|l instead, one finds that 
the pressure even becomes negative above T > 300 MeV. 
Another problem is that for T < Tc their pressure is 
equal to that of a massless gas. However, this is incor- 
rect since the sigma meson is massive and included at 
NLO. Hence one expects a deviation from the ideal-gas 
pressure at T < Tc. Finally, at high temperatures we ex- 
pect that the NLO pressure will become approximately 
equal to the LO pressure because chiral symmetry will be 
restored. This is not the case for the pressure of Ref. 

We briefly comment on the paper by Jain In 
that paper the author is calculating the thermodynamic 



potential to NLO in the 0{N) linear sigma model using 
a high-temperature expansion. Since this approximation 
breaks down at low temperatures, we will refrain from 
comparing with our results. 



CHOICE OF PARAMETERS 



In the preceding sections we have shown plots for 
particular choices of the parameters, namely, A(/i = 
100 MeV) = 30, = 47 MeV and if there is explicit 
symmetry breaking, we take H = (104 MeV)"^. In this 
section we will motivate these choices. For simplicity we 
partly use leading-order calculations for fixing the pa- 
rameters. 

We start with choosing the values for /^r and m-^ to 
be roughly equal to their measured values: = 47 MeV 
(note that our /^r differs from the more conventional defi- 
nition by a factor of 1/2) and m^^ = 138 MeV (the average 
of the measured masses of the 7r°,7r+ and tt"). We will 
use this for choosing our parameter H as follows. Given 
a choice of A at a given scale /i we solve the LO renor- 
malized gap equations H28() and H29() . for a and m^, such 
that = (which is the correct identification at LO). 
For our choice of X(p ~ 100 MeV) = 30, this results in 



(68) H = (104 MeV)-^ 



The choice of A is motivated by considerations on 
the maximal value of the cutoff and the sigma mass. To 
obtain this mass, one has to find the poles of the prop- 
agators in Minkowski space. The physical mass TOph is 
often defined by the solution to the equation 



-I- + ReI](po = imph + e, p — 0, m) — , (70) 



where S is the self-energy. Using Eq. (|19|l and choosing 
/i = mcr, we find that at T = and for H = 



327r2/2 



2 327r' 



A(m<,) l+327r^/A(m„) 



(71) 



Equation (|71(l implies that m^ < VlGTrf-^ ~ 333 MeV, 
which is lower than 600 MeV. For H a, similar bound 
applies. In that case we find that the maximal value of 
the sigma mass can be found by solving the following 
equation for m„ 



where 



A{x) = 



\/l + A^(jn^Jrn^ m^ 



167r/2 1 



(72) 



mi 



TT 



TT I 1+ Vl-4/x 



(73) 
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We find by solving this equation that the maximal value 
of TTicr is equal to 433 MeV. This is also smaller than the 
average measured value of 600 MeV. The reason that we 
find an unphysical bound for the sigma mass could be 
that wc consider Nf = 2 and may miss out on essential 
three-flavor physics. 

The sigma mass turns out to be maximal if A(/i — 
100 MeV) = 80. The problem with this choice of A is that 
in that case the maximal value of the cutoff is 720 MeV. 
This is very low and allows us only to do calculations up 
to around T ~ 50 MeV. Therefore we choose a lower 
value: A(/i = 100 McV) = 30. Using that parameter 
choice Aniax = 19 GcV and the sigma mass is equal to 
256 MeV and 350 MeV in the case of iJ = and iJ = 
(104 MeV)^ respectively. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we have considered the thermodynam- 
ics of the 0{N) linear and nonlinear sigma models to 
NLO in the 1/iV expansion. 

At NLO we have shown that one can renormalize the 
effective potential in a temperature-independent manner 
only at the minimum of the effective potential. This is 
another example of the ambiguity in the definition of off- 
shell Green's functions. A perhaps more familiar exam- 
ple comes from the calculation of the zero-temperature 
effective potential in gauge theories. In this case, the ef- 
fective potential depends on the gauge-fixing condition 
except at the minimum |2^ |25| . By renormalizing 
the NLO effective potential in the minimum we found 



the beta function for A to NLO. This beta function is 
consistent with the perturbative calculation. 

We calculated numerically the pressure for the linear 
and nonlinear sigma model to NLO as a function of tem- 
perature. Our calculations show that for the calculation 
of the pressure 1/iV is a good expansion, even if TV = 4. 
With a relatively realistic choice of the parameters we 
made a prediction for the pressure of QCD at low temper- 
atures. Our results for the pressure disagree significantly 
with the calculations of those in Ref This is due 
to the fact that we are not neglecting zero-temperature 
contributions and that we treat the NLO contribution 
without resorting to any high-energy approximation. 

We also found that in the linear sigma model the 
sigma mass has an upper bound. This bound depends 
only on the parameters and m,r- For a realistic choice 
of these parameters, this implies that the sigma mass is 
smaller than 433 MeV. This does not necessarily have 
consequences for the real sigma meson, since we did not 
take into account the full three-flavor physics. 

Having solved the 0{N) model for the thermody- 
namics, it is natural to apply it to other quantities 
such as spectral functions. The methods developed here 
should also be useful for more complicated models in- 
corporating additional features of low-energy QCD, e.g. 

U{3)a X C/(3)y. 
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